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Abstract. There are 13 equivalence classes of 2D second order quantum and classical 
superintegrable systems with nontrivial potential, each associated with a quadratic algebra 
of hidden symmetries. We study the finite and infinite irreducible representations of the 
, quantum quadratic algebras though the construction of models in which the symmetries 

act on spaces of functions of a single complex variable via either differential operators or 
C"| , difference operators. In another paper we have already carried out parts of this analysis for 

the generic nondegenerate superintegrable system on the complex 2-sphcrc. Here we carry 
it out for a degenerate superintegrable system on the 2-sphere. We point out the connection 
between our results and a position dependent mass Hamiltonian studied by Quesne. We 
also show how to derive simple models of the classical quadratic algebras for superintegrable 
systems and then obtain the quantum models from the classical models, even though the 
, classical and quantum quadratic algebras arc distinct. 
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1 Introduction 



A classical (or quantum) mth order superintegrable system is an integrable n-dimensional Hamil- 



tonian system with potential that admits 2n — 1 functionally independent constants of the 



motion, the maximum possible, and such that the constants of the motion are polynomial of 
at most order m in the momenta. Such systems are of special significance in mathematical 
physics because the trajectories of the classical motions can be determined by algebraic means 
alone, whereas the quantum eigenvalues for the energy and the other symmetry operators can 
also be determined by algebraic methods. In contrast to merely integrable systems, they can 
be solved in multiple ways. The best known (and historically most important) examples are 
the classical Kepler system and the quantum Coulomb (hydrogen atom) system, as well as the 
isotropic oscillator. For these examples m = 2 and the most complete classification and struc- 
ture results are known for the second order case. There is an extensive literature on the subject 
[IllIliailElElEllHligQI^ with a recent new burst 

of activity [SliaESEgEaEElEg 

All such systems have been classified for real and complex Riemannian spaces with n = 2 and 



*This paper is a contribution to the Proceedings of the Seventh International Conference "Symmetry in 
Nonlinear Mathematical Physics" (June 24-30, 2007, Kyiv, Ukraine). The full collection is available at 
|http: / /www.emis.de/journals / SIGMA /symmetry2007.html | 
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their associated quadratic algebras of symmetries computed [15, 29. 301 124^ [31] . For nonconstant 
potentials there are 13 equivalence classes of such stems (under the Stackel transform between 
manifolds), 7 with nondegenerate (3-parameter) potentials and 6 with degenerate (1-parameter) 
potentials [HI [30] . The constants of the motion for each system generate a quadratic algebra 
that closes at order 6 in the nondegenerate case and at order 4 in the degenerate case. 

The representation theory of such algebras is of great interest because it is this quadratic al- 
gebra "hidden symmetry" that accounts for the degeneracies of the energy levels of the quantum 
systems and the ability to compute all associated spectra of such systems by algebraic means 
alone. In principle, all of these quadratic algebras can be obtained from the quadratic algebra of 
a single generic 3-parameter potential on the complex two-sphere by prescribed limit operations 
and through Stackel transforms. However, these limiting operations are not yet sufficiently 
understood. Each equivalence class has special properties, and each of the 13 cases is worthy 
of study in its own right. A powerful technique for carrying out this study is the use of "one 
variable models". In the quantum case these are realizations of the quadratic algebra (on an 
energy eigenspace) in terms of differential or difference operators acting on a space of functions 
of a single complex variable, and for which the energy eigenvalue is constant. Each model is 
adapted to the spectral decomposition of one of the symmetry operators, in particular, one 
that is associated with variable separation in the original quantum system. The possible irre- 
ducible representations can be constructed on these spaces with function space inner or bilinear 
products (as appropriate) and intertwining operators to map the representation space to the 
solution space of the associated quantum system. (There have been several elegant treatments 
of the representation theory of some quadratic algebras, e.g. [Ml [T71 [THl EJ El [Ml [37]. However 
these have almost always been restricted to finite dimensional and unitary representations and 
the question of determining all one variable models has not been addressed.) In the classical 
case these are realizations of the quadratic algebra (restricted to a constant energy surface) by 
functions of a single pair of canonical conjugate variables. 

In [39] we have already carried out parts of this analysis for the generic nondegenerate super- 
integrable system on the complex 2-sphere. There the potential was V = a\/s\ + 02/ 1 s\ + 03/ s| 
where s\ + s\ + s\ = 1, and the one variable quantum model was expressed in terms of difference 
operators. It gave exactly the algebra that describes the Wilson and Racah polynomials in their 
full generality. In this paper we treat a superintegrable case with a degenerate potential. Our 
example is again on the complex 2-sphere, but now the potential is V = ot/s\. Though this 
potential is a restriction of the generic potential, the degenerate case admits a Killing vector so 
the quadratic algebra structure changes dramatically. The associated quadratic algebra closes 
at level 4 and has a richer representation theory than the nondegenerate case. Now we find one 
variable models for an irreducible representation expressed as either difference or differential 
operators, or sometimes both. We show that this system can occur in unobvious ways, such as 
in a position dependent mass Hamiltonian recently introduced by Quesne |37] . 

The second part of the paper concerns models of classical quadratic algebras. Here we inau- 
gurate this study, in particular its relationship to quantum models of superintegrable systems. 
We first describe how these classical models arise out of standard Hamilton-Jacobi theory. 
In [441 [27] we have shown that for second order superintegrable systems in two dimensions 
there is a 1-1 relationship between classical quadratic algebras and quantum quadratic algebras, 
even though these algebras are not isomorphic. In this sense the quantum quadratic algebra, 
the spectral theory for its irreducible representations and its possible one variable models are 
already uniquely determined by the classical system. We make this concrete by showing ex- 
plicitly how the possible classical models of the classical superintegrable system with potential 

V = ct/s| lead directly to the possible one variable differential or difference operator models for 
the quantum quadratic algebra. Then we repeat this analysis for the nondegenerate potential 

V = a\/s\ + 02/^2 + a 3/ s i where the quantum model is essentially the Racah algebra QR(3) 
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and its infinite dimensional extension to describe the Wilson polynomials. Our results show 
that the Wilson polynomial structure is already imbedded in the classical system with potential 
V = a\/s\ + g^/s! + 03/^3, even though this potential admits no Lie symmetries. Thus the 
properties of the Wilson polynomials in their full generality could have been derived directly 
from classical mechanics! 

This work is part of a long term project to study the structure and representation theory 
for quadratic algebras associated with superintegrable systems in n dimensions |23[ [24"1 [25| [26"1 
l27t 135] [36] . The analysis for n = 3 dimensions will be much more challenging, but also a good 
indication of behavior for general n. 



2 The structure equations for S3 

Up to a Stackel transform, every 2D second order superintegrable system with nonconstant po- 
tential is equivalent to one of 13 systems [44] . There is a representative from each equivalence 
class on either the complex 2-sphere or complex Euclidean space. In several papers, in particu- 
lar [15], we have classified all of the constant curvature superintegrable systems, and this paper 
focuses on two systems contained in that list: S9 and S3. The quadratic algebra of the generic 
nondegenerate system S9 was already treated in [39] and we will return to it again in this paper. 
First we study the quadratic algebra representation theory for the degenerate potential S3. This 
one-parameter potential 2-sphere system corresponds to the potential 

V = - 

s 2 ' 

where s\ + s 2 , + s| = 1 is the imbedding of the sphere in Euclidean space. The quantum 
degenerate superintegrable system is 

H = J 2 + Jf + Jf + V(x,y) =H + V, 

where J3 = sid S2 — S2<9 Sl and J 2 , J3 are obtained by cyclic permutations of the indices 1,2,3. 
The basis symmetries are 



L a = Jf + ^, L 2 = i(Ja J 2 + J 2 J,) - 

s 3 I s 3 



where J3 = S2d Sl — sid S2 plus cyclic permutations. They generate a quadratic algebra that 
closes at order 4. The quadratic algebra relations are [H, X] = [H, Lj] = and 



[La, L 2 ] = -(LtX + XLa) - ( \ + 2a ) X. 



[L 1 ,X] = 2L 2 , [L 2 ,X] = -X 2 + H -a, (1) 

'I 

The Casimir relation is 

C = i (X 2 Li + XL X X + L X X 2 ) +L\ + Ll- HL\ + (a + X 2 - 
( 2\ 5a 

+ { a -- 3 ) Ll -ir = - (2) 

We know that the quantum Schrodinger equation separates in spherical coordinates, and that 
corresponding to a fixed energy eigenvalue H , the eigenvalues of X take the linear form 

A n = An + B, 
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where n is an integer, so we will look for irreducible representations of the quadratic algebra such 
that the representation space has a basis of eigenvectors f n with corresponding eigenvalues X n . 
(Indeed, from the analysis of [16] or of [6] the structure equations imply that the spectrum of X 
must be of this form.) We will use the abstract structure equations to list the corresponding 
representations and compute the action of L\ and L 2 on an X basis. Thus, we assume that 
there is a basis {/«}, for the representation space such that 

Xf n = A„/„, Lif n = C(j, n)fj, L 2 f n = ^2 D ti> n )^'- 

3 3 

Here, A, B are not yet fixed. We do not impose any inner product space structure. 

From these assumptions we can compute the action of L\ and L 2 on the basis. Indeed, 

[Li,L 2 ]f n = (C(j, k)D(k, n) - D(j, k)C(k, n)) f„ (3) 

3,k 

[L u X\f n = £>n - Xj)C(j, n)f s , [L 2 ,X]f n = ]T(A n - Xj)D(j, n)f r (4) 
3 3 

On the other hand, from the equations (HD we have 

[L 1 ,X]f n = 2j2D(j,n)f j , (5) 
3 

[L 2 ,X]f n = -2 C(j, n)fj + (-X 2 n + H- a)f n , (6) 
3 

[Li,L 2 ]f n = - ^(A n + \j)Cti,n)fj - Q + 2a\ X n f n . (7) 

3 

Now we equate equations (jlj) with ([5]) or (|6|). For j = n, equating coefficients of f n in the 
resulting identities yields the conditions 

D{n, n) = 0, C(n, n) = ~ A " + g ~ a . 

Similarly, equating coefficients of fj in the case j 7^ n yields 

A{n - j)D(j, n) = -2C(j, n), A(n - j)C(j, n) = 2D(j, n), 

or 

{A 2 (n-j) 2 +A)C(j,n) =0, j + n. 

Thus, either C(j, n and D(j, n) vanish or A 2 (n— j) 2 = —4. We can scale A such that the smallest 
nonzero jump is for j = n ± 1, in which case A = ±2i. By replacing n by — n if necessary, we 
can assume A = 2%. (We also set B = ifi.) Thus the only possible nonzero values of C(j,n), 
D(j, n) are for j = n, n ± 1 and there are the relations 

D(n + 1, n) = —iC(n + 1, n), D(n — 1, n) = iC(n — 1, n). 

Comparing ([3|) and ([7]) and equating coefficients of f n ±2, fn±i, respectively, on both sides of 
the resulting identities, we do not obtain new conditions. However, equating coefficients of f n 
results in the condition 

F n+1 -F n = -(2n + (i) [An 2 + Aan + n 2 + H + a + \ 
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where F n = C(n,n— l)C(n — 1, n). The general solution of this difference equation is 
F„ =n 4 + (2/j - 2)n 3 + Q M 2 - 3/1 + y + | + |) n 2 



1 2 2 



H a 5 

y + 2 + 4 



a r 

Li I n + K, 

P 2 2 4' 



where k is an arbitrary constant. 

To determine k we substitute these results into the Casimir equation (J5]) and set equal to 
zero the coefficients of fj in the expression Cf n = 0. For j / nwe get nothing new. However, 
j = n we find 

1 4 1 , + a + 2 l/2 + # + a 
k = — it u 6 H '—uT a 

1 „ a .ff 2 Ha a 2 

+ -H-\ 1 1 . (8) 

8 8 16 8 16 V ' 

Thus, F n = C(n,n — l)C(n — l,n) is an explicit 4th order polynomial in n. By factoring this 
polynomial in various ways, and re-normalizing the basis vectors f n appropriately via f n — > 
c(n)/ n , we can achieve a realization of the action of L\ and L2 such that 

Lifn = C{n + l,n)/ n+ i + C(n,n)/„ + C(n- l,n)/ n _i, 

L 2 f n = D(n + 1, n)f n+l + D(n, n)f n + D(n - 1, n)/„_i (9) 
and all of the coefficients are polynomials in n. The 4 roots of F n are 



1-H 1 



±-y]2- 4(iJ + q) ± 2 v / l-4(# + a) + 16#a, 



2 4 

so a convenient factorization is 

C{n,n-l) = (n + ^^\ - ^ (l - 2(H + a) + ^1 - A{H + a) + 16Ha 
C(n-l,n) = (n + ^-^-\ - ^ (l - 2(iJ + a) - y/l - A{H + a) + lQHa 

From these expressions and from 

£>(n,n)=0, C(n,n) = i 2n + ^ + H ~ a ^ D(n ± l,n) = ^iC(n ± l,n) 

we see that we can find C, D coefficients in which the the dependence on n is always as 
a polynomial. 

There are raising and lowering operators 

= Li + iL 2 + ^(X 2 - H + a), A = Li - iL 2 + \{X 2 - H + a). 

Indeed, 

A^f n = 2C(n + 1, n)f n+1 , Af n = 2C(n - 1, n)f n -i, 
and [A, A']f n = 2(F n+ i — F n )f n , so [A, A^] is a third order polynomial in X. 
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To get a one-variable model of the quadratic algebra in terms of second order differential 
operators, we can simply make the choices f n (t) = t n , X = H^tA + fi) and define L\ from 
expressions ([9]) via the prescription 

L./.W = [to ( t | + 1,4) + c (4,4) + r'c (4 - 1,4)) /.(.,, (io) 

with a similar procedure for L2. 

In general the irreducible representations that we have defined are infinite dimensional and 
the basis vectors f n occur for all positive and negative integers n. We can obtain representations 
bounded below, and with lowest weight [i for — iX and corresponding lowest weight vector /o, 
simply by requiring Fq = 0, which amounts to setting k = 0. For convenience we set a = 1/4— a 2 . 
Then we have 

F n = C(n,n- l)C(n - l,n) = n(n + (i - l)(n + + a)(n - a). (11) 
Since k = 0, (|8|). must be a solution of this quadratic equation: 

H = -( Ji -l + a) 2 + ± (12) 
A convenient choice is 

C(n — 1, n) = n(n + fx — 1 + a), C(n + 1, n) = (n + fj)(n + 1 — a), 
C(n,n) = 2n 2 + 2n/i — fia + a + fi — — . 



If /i is not a negative integer then this bounded below representation is infinite dimensional. 
However, if there is a highest weight vector f m then we must have F m+ i = 0, or /i = — m, 
m = 0, 1, . . . . Thus the finite dimensional representations are indexed by the nonnegative integer 
m and the eigenvalues of — iX are m, m — 2, . . . , — m. The dimension of the representation space 
is m + 1. 

At this point it is worth pointing out that all of the finite dimensional, infinite dimensional 
bounded below, and general infinite dimensional irreducible representations and models of the 
quadratic algebra associated with the superintegrable system S3 are of direct interest and ap- 
plicability. A similar argument was made in [39] where we gave examples of various analytic 
function expansions and distinct unitary structures associated with one superintegrable sys- 
tem. The original S3 quantum system is given in terms of complex variables with no specific 
inner product structure imposed. One could use the representation theory results to describe 
eigenfunction expansions simply in terms of analytic functions. If one wants an inner product 
structure or bilinear product structure, it is merely necessary to impose the structure on a single 
eigenspace of H, and there are a variety of ways to do this. For example, one could restrict 
the complex system to the real sphere and impose the standard inner product for that case. 
Alternatively, one could restrict to the real hyperboloid of one sheet, or the real hyperboloid of 
two sheets. In all cases the models of the irreducible representations are relevant, though not 
necessarily to one special case, such as the real sphere with the standard inner product. While 
we have no direct proof that all models of irreducible representations of quadratic algebras ob- 
tained in this way lead to representations for some version of the original superintegrable system, 
we have no counterexamples. For a deeper analysis we need to construct intertwining operators 
that relate basis functions for the model with eigenfunctions of the quantum Hamiltonian. 
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3 Differential operator models 

A convenient realization of the finite dimensional representations by differential operators in 
one complex variable is 

d 2 d 
Li = (t 3 + 2t 2 + t) —J + ((2 - a - m)t 2 + 2(1 - m)t + a - m) — 

Lib Lbu 

+ m(a — l)t + a(m + 1) — m — — , X = i(2t-^ — m), 



d 2 d 
L 2 = i (-t 3 + t) —T +i((a + m- 2)t 2 + a - m) — - im(a - l)t. (13) 
dt z dt 



This model is also correct for infinite dimensional bounded-below representations, except that 
now the lowest weight is fj, = —m where m^0,l,2,... is a complex number. The raising and 
lowering operators for the model are 

= 2t 3 ^r + 2(2 -a- m)t 2 4- + 2m(a - l)f, A = 2t^-r + 2(a - m)4- 
dt z dt dt z dt 



In the finite dimensional case, for example, the eigenvalues of L\ are 

Xn = a 2 - \ - ( n - a + - J , n = 0, 1, . . . ,m, (14) 



4 V 2 . 
and the corresponding unnormalized eigenfunctions are 



:i+t)v*i n ~ a n - m -,-t). 

a — m 



Now, motivated by the quantum mechanical system on the real 2-sphere, we impose a Hilbert 
space structure on the irreducible representations such that L\ and L 2 are self-adjoint and X is 
skew adjoint: 

(Ljfn, fn') = (fn,Ljf n >), j = 1,2; (X f n , f n >) = — (f n , X f n r) . 

Writing cj) n = k n f n where 4> n has norm 1, we have the recursion relation 
2 _ (n - 1 + fi)(n - a) u2 

"n- n{n _ 1 + fI + a) Kn-l- 

For infinite dimensional bounded below representations k 2 must be positive for all integers 
n > 0, and we normalize ko = 1. Thus 

k 2 = Mn(l ~ a)n 
n\(a + n) n 

For finite dimensional representations we have fi = —m. Normalizing k$ = 1, (possible for 
a < 1 or for a > m), we find that an orthonormal basis in the one variable model is given by 
<j) n (t) = k n fn(t) = k n t n , n = 0, 1, . . . , m where 



/ (-m) n .(l - a) n _ / m!(l - a) n (l - a) m -n 
Y n\(-m + a) n \j (1 - a) m n\(m - n)! 

Note the reflection symmetry \\f n \ \ = 1 1 fm-n \ \ ■ 
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To derive a realization of the Hilbert space for the differential operator models of the finite di- 
mensional and infinite dimensional bounded below unitary representations in terms of a function 
space inner product 



(p,q)=K J J p(t)q(t)p(tt)dtdt, 



where p, q are polynomials and K is a normalization constant, we use the formal self- and 
skew-adjoint requirements and obtain a differential equation for the weight function: 

(-C 2 + C)^P + (-/^-a + l + (-l + M-a)C)^ + (-2 + / u-2a + a^(C) = 0, 

where ( = it. The solution that vanishes at ( = 1 for a < 1/2 and is integrable at ( = for 
a + p. > — 1 is 

/j+3a-Q _| -i _/i+3a+Q 



( ^+3a-Q _| -i /i+3a+Cj , -i \ 

2-2a 2 



where Q = -\f a 2 + (2/i — 8)a + /U 2 + 4/i — 8. (Note that the integral is an even function of Q.) 
At ( = this function has a branch point with behavior ( a +A\ We write t = re ld , t = re~ ld , 
( = r 2 and choose our contours of integration for the inner product as the unit circle \e l ®\ = 1, 
i.e., < 9 < 2ir and, in the complex ("-plane, a contour that starts at ( = 1 and travels just 
above the real (-axis to circle ( = once in the counterclockwise direction and returns to ( = 1 
just below the real (-axis. We require that (1, 1) = 1. By choosing a regime where a + p > — 1 
we can shrink the (-contour about ( = so that the norm takes the form 

(1, 1) = -4TrK 1 e in{a+ll) sm[-ir(a + p)] [ pi{£)<K 

Jo 



= AttK, sinWa + a)} r ( 2 ~ 2a ) r ( a + ^ + ^ 

where we have integrated term-by-term and then made use of Gauss' Theorem for the summation 
of 2-Fi(l)- This gives us the value for K\ such that (1, 1) = 1. Now, the result extends for the 
original contour by analytic continuation. This defines a pre-Hilbert space inner product that 
then can be extended to obtain a true Hilbert space. 

The contour integral for the inner product obtained in the previous paragraph requires 
Re a < 1 for convergence, and this doesn't hold for some of the unitary irreducible representa- 
tions defined above. Accordingly, we consider a second solution of the weight function equation. 
The solution that vanishes at ( = for a + p, > and is integrable at ( = 1 for a < 1 is 

( /j+3a— Q fi+3a+Q 
\i + a + 1 

At ( = 1 this function has a branch point with behavior (1 — () 1_2a . We write t = re 10 , t = re~ ld , 
Q = r 2 and choose our contours of integration for the inner product as the unit circle \e ld \ = 1 
and, in the complex (-plane, a contour that starts at ( = and travels just below the real (-axis 
to circle ( = 1 once in the counterclockwise direction and returns to ( = just above the real 
(-axis. This integral converges for Re (a + fi) > —1. We require that (1, 1) = 1. By choosing 
a regime where a < 1 we can shrink the (-contour about ( = 1 so that the norm takes the form 



(1, 1) = 47r/f 2 e^ (2a - 1) sin[7r(2a - 1)] [ 

Jo 



P2(()d( 
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K / M+3a+Q ti+3a-Q \ 



/U + a + 1 

-47ris:2e i7r ^ a - ij sin[7r(2a - 1 



i*(2«-l) ^_ r _,o„ i m r ( 2 " 2 «) r ( a + M + 1) 



'r(2 - a ~^+Q )r(2 - a ~^~° )" 



This gives us the value for .fT such that (1,1) = 1, and the result extends by analytic continuation 
to all values of a, \x for which the original contour integral converges. 

Thus we have an explicit pre-Hilbert function space inner product for each of our differential 
operator models. In the finite dimensional case we have the reproducing kernel function 



n=0 

In the infinite dimensional bounded below case we have the reproducing kernel function 



*(«,*) = f>(t)0„(*) = 2^ /. n 1 a ;^) 



n=0 

which converges as an analytic function and in the Hilbert space norm for \s\ < 1. Here, 
In each case (f(t),6(t,s)} = f(s) for / in the Hilbert space. 



4 Difference operator models 

There are also difference operator models for the representations of the S3 quadratic algebra. 
We first give the details for the finite dimensional representations indexed by the nonnegative 
integer m. Here the operator L\ is diagonalized: 

Li = -A(i) + a - -, X(t) =t{t-2a + l), 

_. x = (t-2a + l)(t-m) Tl _ t(t + m-2a + l) T _ 1 
2t-2a + l 2t-2a + l 

(t-a + l)(t-2a + l)(t-m) rrl , t{t - a){t + m - 2a + 1) _ x 

ho = 1 H 1 , 

2t-2a + l 2t-2a + l 

where T k is the difference operator T k f(t) = f(t + k). The basis functions are f n {t) = 
(—l) n p n (X) where 

Pn{X{t)) = 3 F 2 ( 1 

y — m I — a 

Here f n is a polynomial of order n in the variable A(t), a special case of the family of dual Hahn 
polynomials \45\ p. 346]. These polynomials are orthogonal with respect to a measure with 
support at the values t = 0, 1, . . . , m, in agreement with equation (|14p for the eigenvalues of L\. 
Indeed, we have (for a < 1) 

^ (l-2a) t (3/2-a),(-m-l) t (-l) f (2 - 2a) m (a - m) n n\ 

^ (1/2 - a) t {2 + m - 2a) t t\ (1 - a) m (l - a) n (-m) n 
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For the infinite dimensional, bounded below, case we have 
1 



L 1 = t 2 + a 2 



4' 



_ iX = (1/2- a- it)(u + a- 1/2 -it) ^ _ (1/2 - a + it)(a + a - 1/2 + it)) T _ 



2t 

.(1 - 2it)(l/2 -a- it)(fi + a -1/2- it) i 
4t 

.(l + 2ii)(l/2 - a + it)(u + a- 1/2 + it)) 

- 1 1 

4t 



2t 



(15) 



The basis functions are f n (t) = (—l) n s n (t 2 ) where 



Sn{t 2 



3F2 



—n t; — a + it i — a — it 
H 1-a ' ' ' 



Here f n is a polynomial of order n in the variable t 2 , a special case of the family of continuous 
dual Hahn polynomials [45^ p. 331]. The orthogonality and normalization are given by 



1 r 

2^ Jo 



I\l/2 - a + it)T(ji + a - 1/2 + it)T(l/2 + it) 



T(2it) 

r(n + /i)r(n + 1 — a)T(n + fx + a)n! 

(^I(l-«)n| 2 



S n (t 2 ) Sn r(t 2 )dt 



where fi > 1/2 — a > 0. 

In summary, we have found the following possibilities for bounded below irreducible repre- 
sentations such that L±, L2 are self-adjoint and X is skew adjoint, together with associated one 
variable models. (Here, uq is a positive integer.) 



representation 


parameter range 


model 


finite dimensional 
inf. dim. bdd. below 
inf, dim. bdd. below 
inf. dim. bdd. below 


fi = —m, m = 0, 1, 2, . . . 
either a<lora + /i>0 

A* > 
a < 1 and a + fx > 

> fi = -n + t, t e (0, 1) 
a = no + s, s G (0, 1) 

> // = -n + t, t e (0, 1) 
— f < a < 1 — t 


differential operators 
difference operators 

differential operators 
difference operators 

differential operators 
differential operators 



5 Quesne's position dependent mass (PDM) system 
in a two-dimensional semi-infinite layer 

In [37] Quesne considered a superintegrable exactly solvable position dependent mass (PDM) 
system in a two-dimensional semi-infinite layer. Her system is equivalent via a gauge transfor- 
mation to a standard quantum mechanical problem on the real 2-sphere with potential of the 
form S3. Indeed, in Quesne's paper we are given the Hamiltonian 

—Hn = cosh 2 qx(d 2 + d 2 ) + 2q cosh qx sinh qxd x + q 2 cosh 2 qx — - — ^—r . 

sinh qx 
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We adopt coordinates on the unit sphere as 

sin qy cos qy 

si = , s 2 = : , s 3 = tanhgx, 

cosh qx cosh qx 

where s\ + s 2 + s| = 1 and the metric is ds 2 = q 2 (dx 2 + dy 2 )/ cosh 2 qx. The Laplacian becomes 

cosh 2 <?x 2 2 
As = ^ + 

In these coordinates, the degenerate superintegrable system S3 becomes 

_ cosh 2 qx 2 2 \-a 2 

Hs -^^ {d * + d y ) + ^Wq-x 

By a gauge transform Hq = (cosh qx)~ 1 Hg cosh gx, we get 

if 1 - a 2 1 1 

iio = -5- \ cosh 2 qx(d 2 + d 2 ) + 2g cos gx sinhgx<9 x + q 2 cosh 2 gx + q 2 — — 5 > + — — a 2 . 

q I y sinh qx I 4 

Thus we have i?Q = —q 2 Ho + (7 2 (l/4 — a 2 ) = — g 2 (cosh qx)~ 1 Hs cosh gx + q 2 (l/A — a 2 ), with 
1/4 — a 2 = —k(k — 1) which has solutions k = a + 1/2 or k = —a + 1/2. Since k is assumed 
positive and a is required to be less than 1, we take a < 0. 

Suppose we find an eigenvector for Hs with eigenvalue As, call it v\ s . Then Aq will be the 
eigenvalue of v\ Q for Hq. We have the transformations v\ Q = v\s/cosh.qx and Aq = —q 2 Xs + 
g 2 (l/4— a 2 ) = —q 2 Xs—q 2 k(k—l). Checking the two eigenvalues, we have As = — (fi— l+a) 2 + l/4 
and Aq = q 2 (N + 2)(N + 2A: + 1) . We note that these two values coincide when — \i = m = N + 1 
with m an integer. 

Using the above calculations and the eigenfunctions given in the paper, we can obtain eigen- 
functions for the S3 case as 

^ As = iV^ ) (tanhgx)' a+ ^(coshgx)^- 1 p- a ^ +1 (-tanh 2 gx)x£(y), 
or in coordinates on the sphere 

= N { n %3)- a+ ksi + s 2 )^P-^(-s 2 ) Xe (y), 

where m = 2n + £+l, and Xiiu) = sin[(^+ or cos[(£+l)qy]. We can rewrite these by noting 
l/cosh 2 gx = 1 — tanh 2 qx so that we can write sin qy = s\/ 'yl — s 2 and cos qy = s<ij \j 1 — s 2 , 
then we obtain 

Xe(y) = a n T e+ i 82 + b n Sl U e 



V4+4J V4+4 \V4+4J 

where and Up are the Chebyshev polynomials of the first and second kind, respectively. 

Quesne found the S3 quadratic algebra (which closes at order 4) but did not use it for 
spectral analysis purposes because her problem involved a boundary condition that broke the 
full quadratic algebra symmetry. Instead she considered her system as a special case of S9 
and used the more complicated S9 symmetry algebra that closes at order 6 to find the finite 
dimensional representations. (Note that although the 1-parameter S3 potential is a limit of the 3- 
parameter S9 potential as two of the parameters go to 0, a discontinuity occurs in the structure 
of the quadratic algebra. A first order symmetry appears and the number of second order 
symmetries jumps from 3 to 4.) Quesne's point of view has merit, but it complicates the spectral 
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analysis of the problem, since the only one- variable model is in terms of difference operators and 
Racah polynomials. From our vantage point of one one variable differential operator analysis 
for the model, Quesne's boundary conditions amount to decomposing an irreducible subspace 
corresponding to an m-dimensional representation into a direct sum of even and odd parity 
subspaces V + , V~ . (Indeed her boundary conditions require choice of Xiiv) i n the cosine 
form for I even and in the sine form for i odd.) Let P be the operator Pf(t) = t m f(l/t). 
Since P 2 = I and ||/ n || = ||/ m - n ||, it is clear that P is unitary. We define unit vectors 
= 2-V2(^ + (-i)^^-) and = 2~ 1 /2(^ _ (_i)-^) for i = 0, 1, ... , [m/2]. Then 

for m = 2k the vectors I = 0, . . . , k form an orthonormal basis for and the vectors , 
I = 0, . . . , k — 1 form an on basis for V^. For m = 2k — 1, the vectors I = 0, . . . , k — 1 form 
an ON basis for and the vectors <3?^~, £ = 0, . . . , k — 1 form an orthonormal basis for V~. 
These basis vectors are very easily expressible in terms of the one variable differential operator 
model, where they are sums of two monomials. The basis used by Quesne corresponds to the V~ 
subspaces. Thus our models can be used to carry out the spectral analysis for this PDM system, 
and they yield a simplification. 

5.1 Classical models for S3 

Now we describe how the methods of classical mechanics lead directly to the quantum models. 
The classical system S3 on the 2-sphere is determined by the Hamiltonian 

n, ^2,^2,^2, a ( s i + s 2 + s i) 
ri = J\ + J2 + H 2 ' 

where J\ = S2P3 — S3P2 and ^2, Jj, are cyclic permutations of this expression. For computational 
convenience we have imbedded the 2-sphere in Euclidean 3-space. Thus we use the Poisson 
bracket 

3 

{F,G} = Y,(-ds^d Pi g + d Pi Td Si g) 
1=1 

for our computations, but at the end we restrict to the sphere sf + s| + s| = 1. The classical 
basis for the constants of the motion is 

s 2 s s 

C\ = J 2 + a-|, £2 = J1J2 — ol — X = Jz- 

The structure relations are 

{X,d} = -2C 2 , {X,C 2 } = 2C 1 -H + X 2 + a, {£1, £ 2 } = -2(A + a)X, (16) 
and the Casimir relation is 

C\ + L\ - CiH + C1X 2 + aX 2 + ad = 0. (17) 

From the results of [56] we know that additive separation of variables in the Hamilton-Jacobi 
equation TL = E is possible in subgroup type coordinates in which X, C\ or S = 2{C\ — 1L2) — 
TL + X 2 , respectively, are constants of separation. This corresponds to two choices of spherical 
coordinates and one of horospherical coordinates, respectively. We seek two variable models for 
the Poisson bracket relations (|16p . (|17p . There is also separation in ellipsoidal coordinates (i.e., 
non-subgroup type coordinates) but we will not make use of this here. 

The justification for these models comes from Hamilton-Jacobi theory. The phase space for 
our problem is 4-dimensional. Thus it is possible to find canonical variables TL, X, Q, V such 
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that {2, TC} = {V, Q} = 1 and all other Poisson brackets vanish. In terms of TC and the other 
canonical variables the Poisson bracket can be expressed as 

G) = -d n Td x Q + dxTduQ - d Q Td v Q + d v TdgQ. (18) 

(As follows from standard theory |47j one can construct a set of such canonical variables from 
a complete integral of the Hamilton-Jacobi equation. Our 2D second order superintegrablc 
systems are always multiseparable, and each separable solution of the Hamilton-Jacobi equation 
provides a complete integral. Thus we can find these canonical variables in several distinct ways.) 
Now we restrict our attention to the algebra of constants of the motion. This algebra is generated 
by TC, £i, C±, X, subject to the relation (fT7|) . Thus, considered as functions of the canonical 
variables, the constants of the motion are independent of TC. If we further restrict the system 
to the constant energy space TC = E then we can consider TC as non varying and every constant 
of the motion T can be expressed in the form J-(E, Q, V). This means that the Poisson bracket 
of two constants of the motion, J 7 , Q can be computed as 

{JF, g] = -d Q Td v Q + dvTdgQ. 

Thus all functions depend on only two canonically conjugate variables Q, V and the parameter E. 
This shows the existence and the form of two variable models of conjugate variables. However 
the proof is not constructive and, furthermore, it is not unique. Two obtain constructive results 
we will use the strategy of setting Q equal to one of the constants of the motion that corresponds 
to separation of variables in some coordinate system, and then use (I18p for the Poisson bracket 
and require that relations (|16|) . (|17|) hold. In order to make clear that we are computing on the 
constant energy hypersurface expressed in canonical variables we will use a different notation. 
We will set Q e = c,Ve = P so, F(H, Q,V) = /(c,/3), G(TC, Q,V) = g(c,fi), and 

{F, G}e = {/, g} = -djdpg + dpfd c g. 

For our first model we require X = Xe = c. Substituting this requirement and TC = E into 
the structure equations we obtain the result 

1,„ O v 1 



L x = -(E - c 2 - a) + - ^c 4 - 2c 2 (E + a) + (E - a) 2 sin 2/3, (19) 



X = c, L 2 = -^c 4 - 2c 2 (E + a) + (E - a) 2 cos 2^. 

In this model, and in all other classical models, (3 is not uniquely determined: we can replace it 
by = j3 + k(c) for any function k{c) and the variables c and j3' remain canonically conjugate. 
For a second model we require L\ = {C\)e = c and proceed in a similar fashion. The result is 

II: Li = c, L 2 = y/c{E -c-a) sin(2Vc + a(3), (20) 

X = ^J C ( E ~ + C ~ a) "cos(2V^+^0. 

For the third and last model we need to diagonalize the symmetry S = 2{C\ — iL-i) — TC + X 2 
corresponding to separation in horospherical coordinates. For this it is convenient to rewrite the 
structure equations ([To]) . (fT7|) in terms of the new basis S, L\ + i£,2, X- 

{cS, X, } = 2i(S + a), {S, d + iC 2 } = -2iX(S - 2X 2 + 2TC + 3a), 
{£i + iC 2 , X] = -i (X 2 + 2(£i + iC 2 ) -H + a). 

The Casimir relation is 



-2S(d + iC 2 ) - SX 2 + X 4 + TCS- 2X 2 TC + TC 2 -a (2(£i + iC 2 ) + S + 3X 2 + TC) = 0. 
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For model III we set S = c and obtain 



III: S = c, X = -2i(c + a)P, 



Li + iL 2 = 8(c + a) 3 /3 4 + 2(c + a)(3a + c + 2£)/? 2 - ( ° + ^ 



2(c + a) 



5.2 Classical model — > quantum model 



What have we achieved with these classical models? For one thing they show us how to paramete- 
rize the constants of the motion and exhibit their functional dependence. More important for 
our purposes, they give us a rational means to derive the possible one-variable quantum mo- 
dels. This may seem surprising. How can classical mechanics determine quantum mechanics 
uniquely? How can structures such as the Wilson family of orthogonal polynomials, containing 
the Hahn polynomials, be derived directly from classical mechanics? The point is that the struc- 
tures we are studying are second order superintegrable systems in 2D. In papers [30], HU [27] it 
has been shown that there is a 1-1 relationship between the quantum and classical versions for 
such systems, for all 2D Riemannian spaces. (Similarly there is a 1-1 relationship in 3D for 
nondegenerate potentials on conformally flat spaces.) The structures are not identical, since as 
we can see from the examples in this paper, the structure relations in the classical and quantum 
cases are not identical; there are quantum modifications of the classical equations. Although 
we know of no direct prescription for their determination, nonetheless the quantum structure 
equations are uniquely determined by the classical structure equations. Given a second system 
of second order constants of the motion we write down the corresponding quantum system via 
the usual correspondence, where products of classical functions are replaced by symmetrized 
quantum operators, and generate the quadratic algebra by taking repeated commutators. Even 
order classical symmetries correspond to formally self-adjoint quantum symmetries, and odd 
order classical symmetries correspond to formally skew-adjoint quantum symmetries. (This re- 
lationship no longer holds for third order superintegrable systems |21]l40j.) We will demonstrate 
here how to get quantum models from the classical ones that we have derived. 

The basic prescription for the transition from the classical case to the operator case is to 
replace a pair of canonically conjugate variables c, (5 by c — ► t, f3 — > <9j. (There is no obstruction 
to quantization for second order superintegrable systems.) Once an appropriate choice of (3 is 
made in a classical model, we can use this prescription to go to a differential operator model of 
the quantum structure equations. In particular model III above suggests a operator model such 
that S is multiplication by c, X is a first order differential operator in c and L\ + iL 2 is a fourth 
order differential operator. The result, whose existence is implied by the 1-1 classical/quantum 
relationship for second order superintegrable systems, is 

III : S = t, X = -2i(t + a)8t + 2i, (21) 
L x + iL 2 = 8(t + afdf + 2(t + a) (3a + t + 2E + 9)df - 2(t + 5a + AE + 18)d t 
f E a\ E 2 + 2E(9 - a) + (a + 12)(a + 6) 
+ \ +1 2~2J + 2(t + a) ' 

The leading order differential operators terms agree with the classical case but there are lower 
order correction terms needed to correct for the noncommutivity of t and <9j. We can realize 
various irreducible representations of the quadratic algebra by choosing subspaces of functions 
of t on which the operators act. This model agrees with (|10p . (Ilip . fjl2[) in the case where 
C{n — 1, n) = 1 and C(n + l,n) is fourth order. However, there we had a space spanned by 
a countable number of eigenvectors of the skew-adjoint symmetry X whereas here we want the 
spectral decomposition of the self-adjoint symmetry S to govern the model. This forces L 2 to 
be skew-adjoint and X to be self-adjoint. Thus, though the differential operators are formally 
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the same, the Hilbert spaces and the spectral analysis are different. All the representations are 
infinite-dimensional. One class can be realized by closing the dense subspace of C°° functions 
with compact support on < t < oo where the measure is dt/t. The the spectrum of S 
is continuous and runs over the positive real axis. Here X also has continuous real spectra 
covering the full real axis. In particular the generalized eigenfunction of X with real eigenvalue A 
is proportional to and /i is pure imaginary. Thus the spectral analysis of X is given by 
the Mellin transform. There is a similar irreducible representation defined on — oo < t < 0. 
By a canonical transformation we can also get models of these representations in which both 
C{n — 1, n) = 1 and C(n + 1, n) are second order. (We shall illustrate this explicitly for model I.) 
Then the spectral decomposition of S is given by the Hankel transform. Since these particular 
eigenspaces of H admit no discrete spectrum for any of the symmetries of interest, we shall not 
analyze them further. 

Now we consider model I, (I19D . Due to the presence of trigonometric terms in (3 we cannot 
realize this as a finite order differential operator model. However, we can perform a hodograph 
transformation, i.e. use the prescription /3 — > t, c — ► — dt to realize the model. This would seem 
to make no sense due to the appearance of functions of c under the square root sign. However, 
before using the prescription we can make use of the freedom to make a replacement j3' = (3+g(c) 
which preserves canonical variables. We choose 



e -W _> e-^/yV - 2c 2 (E + a) + {E- a) 2 
but leave c unchanged. Then we find 

L x = l(E - c 2 - a) - - [(c 4 - 2c 2 (E + a) + (E - a) 2 ) e 2i/3 - 1 



c i -2c 2 {E + a) + {E-a) 2 )e 2t > j + l , (22) 



with X as before. Now we apply the quantization prescription j3 — > t, c — > —dt and obtain 
a model in which both L\ and L2 are fourth order and X is a first order differential operator. 
This is, in fact, identical to within a coordinate change to model (|2ip . One might also try 
to obtain a difference operator model from (I22p with the replacement c — > t, (3 — > dt, so 
that e 2 *^ would become a difference operator. However, this difference operator quantum model 
is equivalent to what would get from the (3 — > t, c — ► —dt model by taking a Fourier transform. 
Thus we don't regard it as new. 

There is an alternate way to obtain a quantum realization from model I. We use the fact that 

c 4 - 2c 2 (E + a) + (E- a) 2 = (c 2 - (E + a)) 2 - 4a 

and set 



arctan 



V-4a 



-{E + a 
Now we let 2/3 — ► 2/3 + (j) to obtain 

L\ = \{E-c 2 -a) + - ((c 2 - {E + a) 2 ) sin 2/3 + 2iy/acos 2(3) , 

L 2 = - ((c 2 -(E + a) 2 ) cos 2/3 - 2i^sin 2(3) , X = c. 

Now the prescription (3 — ► t, c — > —dt leads to a quantum realization of L\, L2 by second order 
differential operators. Indeed 

L\ = ^(cos(2t) - \)d 2 t - 8i£sm(2t)d t +(~ + 64£ 2 + 8i£ - \ - % J cos(2t) + 1 



2 y v ' ' 1 " v/l V2 " "42/ v ' 2 
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L 2 = i sm{2t)d 2 + 8^cos(2t)d t + (-j + 64£ 2 + 8^ - i - |) cos(2*J + 

Here £ is arbitrary and can be removed via a gauge transformation. The change of variable 
r = e 2tt reduces this model to the form (|13|> . This shows that the flexibility we had in construc- 
ting differential operator models from the abstract representation theory by renormalizing our 
basis vectors f n is replaced in the classical model case by appropriate canonical transformations 
c — ► c, (3 — > P + g(c). In either case there is essentially only one differential operator model that 
can be transformed in various ways. 

It is clear that model II cannot produce finite order differential operator realizations of 
the quantum quadratic algebra, due to the intertwining of square root dependence for c and 
exponential dependence for (3. However, it will produce a difference operator realization via 
Taylor's theorem: e adt f(t) = f(t + a). To show this explicitly we make a coordinate change such 
that 1y/ c + ad c = dc in (|2Q|) . which suggests realizations of the quantum operators in the form 

Lxf(t) = (t 2 - a)f(t), Xf(t) = h(t)f(t + + m(t)f(t - i), 

L 2 f(t) =~(i + 2t)h(t)f(t + %)+ + 2t)m(t)f(t - i). (23) 

A straightforward computation shows that the quantum algebra structure equations are satisfied 
if and only if 

,., . ., 1 (a - t 2 - it)(t 2 +it - E) 

h(t)m(t + i) = -r ~ -7^, 1 ■ (24) 

4 t{t + l) 

Since a = —a 2 + \ and E = — (// — 1 + a) 2 + j for bounded below representations, we can 
factor (|24p simply to obtain 

h(t)m(t + i) = -4^) (t + \ + «a) (t + \ - ia 

t + ^ + ifi + ia^J \ t + ^ — ifj, — ia^J . (25) 

Note that only the product (|24p is determined, not the individual factors. Thus we can 
choose h(t), say, as an arbitrary nonzero function and then determine m(t) from (|24p . All 
these modifications of the factors are accomplished by gauge transformations on the represen- 
tation space: f(t) = p(t)f(t) where p(t) is the gauge function. If we choose the factors in the 
form 

_ , (| - a-it)(p + a - \ -it) _ _ , (| - a + it)(p + a - \ + it) 

~ 1 2t ' rn \ t ') ~ 1 2f ' 

then we we get exactly the model (|15p . The finite dimensional model is related by the simple 
change of variables t — > i(t — a + 1/2), p = —m. In any case, there is only a single solution of 
these equations, up to a gauge transformation. 



5.3 The classical model for S9 



This is the system on the complex sphere, with nondegenerate potential 
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where s\ + s 2 + s i = 1- The classical S9 system has a basis of symmetries 

£l = Ji 2 + o 2 ^| + a3^|, £2 = J 2 2 + a3^ + ai%, Cz = Jl + a l S -% + a 2 S -±, (26) 

where H = £1 + £2 + £3 + ai + a 2 + 03 and the Ji are defined by J% = s\p S2 — S2p Sl and cyclic 
permutation of indices. The classical structure relations are 

{Ci,n} = 8d(H + ai + a 2 + 03) - 8£ 2 - 16£i£ 2 - 16a 2 £ 2 

+ 16a 3 (W + ai + a 2 + a 3 - £1 - £ 2 ), 
{£2, K} = -8£ 2 (H + 01 + a 2 + a 3 ) + 8£^ + 16£i£ 2 + 16ai& 

- 16a 3 (H + ai + a 2 + a 3 - £1 - £ 2 ), 

with {£1, £2} = 1Z and 

K 2 - 16£i£ 2 (H + 01 + a 2 + a 3 ) + 16£ 2 £ 2 + 16£i£^ + 16ai£ 2 

+ 16a 2 £2 + l6a 3 (H + at + a 2 + a 3 ) 2 - 32a 3 (W + a x + a 2 + a 3 )(£i + £ 2 ) 
+ 16a3£ 2 + 32<23£i£2 + 16a3£ 2 — 64aia 2 a 3 = 0. 

Taking C\ = c, H = E with c, /3 as conjugate variables, we find the model 

r 1/ ^ 9 ^ P x («2-Q3)(«i + 2a2 + 2a 3 + £;) 

£2 = 7j(ai + 2a 2 + E - c — — — r 

2 2(c + a 2 + a 3 ) 

V 7 (4aia 2 + 4aia 3 + 2c(£' + ai +a 2 + a 3 ) +4cai — (£ + ai + a 2 + a 3 ) 2 - c 2 )(4a 2 a 3 - c 2 ) 

2(a 2 + a 3 + c) 

x cos(4/3V«2 + a 3 + c). (27) 

This suggests a difference operator realization of the quantum model. 

In the quantum case the symmetry operators L\, L 2 , L 3 are obtained from the corresponding 
classical constants of the motion (|26p through the replacements J k — > J k where the angular 
momentum operators J k are defined by J 3 = x\d X2 — x 2 d xi and cyclic permutation of indices. 
Here H = L\ + L 2 + L 3 + a\ + a 2 + a 3 . The quantum structure relations can be put in the 
symmetric form 

[L t ,R] = 4{Li, L k } - 4{L t ,Lj} - (8 + lGajLj + (8 + 16a k )L k + 8( aj - a k ), 

R 2 = -{L U L 2 , L 3 } + -(16ai + 12) if - (16a 2 + 12)1% - (16a 3 + 12)L| 
6 

52 1 
+ —({L U L 2 } + {L 2 , L 3 } + {L 3 , L x }) + -(16 + 176ai)Li 

1 1 32 

+ -(16 + 176a 2 )L 2 + -(16 + 176a 3 )L 3 + y («i + a 2 + a 3 ) 

+ 48(aia 2 + a 2 a 3 + a 3 ai) + 64aia 2 a 3 . 

Here i, j, A; are chosen such that = 1 where e is the pure skew-symmetric tensor, R = [L x , L 2 ] 
and {L\,Lj} = LiLj + LjLi with an analogous definition of {L\, L 2 , L 3 } as a sum of 6 terms. 
In practice we will substitute L 3 = H — L\ — L 2 — a\ — a 2 — a 3 into these equations. 

Proceeding exactly as in the S3 case (j23[) . (|24p. (|25p . we find that the difference operator 
analogy of (I27p for the quantum quadratic algebra is 



Li = 4i 2 -i + /3 2 + 7 2 , 
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L 2 = h(t)T + m^T- 1 + t(t) 



-4a 2 



8a - 4 + 4£ 2 + m(a + l)t + 16i 2 ](/3 + 1 + 7 - 2it)(/3 - 1 - 7 + 2it) 



+ 



-4a 2 



-2t 2 



1024t(t + i)(2t + i) 2 
4 + 8a + 4£ 2 + 16i(l - a)t + 16i 2 ](/3 + 1 



7 - 2it){d - 1 + 7 + 2it)t L + T~ 



1^2 1 2 1 2 (l 2 - (3 2 )(-4a 2 + 4£ 2 ) 

-3 H — a H — 'Y 4- — — — 

2 P 2 + 2 7 + 8(1 + 4*2) 



where 



° 1 = 4 



a" 



«2 



4 P ' 



fl3 = 4 



H = - 

4 



£ 2 . 



The quadratic terms factor into simple linear terms, and just as in the S3 case, it is only £(t) 
and the product h{t)m(t + i) that is uniquely determined. We can change the individual fac- 
tors h(t), m(t) by a gauge transformation. With the change of variable t = ir and a gauge 
transformation to an operator with maximal symmetry in r, we obtain the standard model 



h{t) = h(r) 

m(t) = m(r) = 
_ £ + a + 1 



(A + t){B + t){C + t)(D + t) 
4t(t + 1/2) 
(A - t)(B - t)(C - t)(D - t) 



D 



4t(t- 1/2) 
E-a + 1 



C 



/3 + 7 + 1 



D 



It follows that L 2 = 1i{t)E +1 +m(T)E- 1 +1{t) is a linear combination of L\ and the difference 
operator whose eigenfunctions are the Wilson polynomials, just as found in [39]. Here E s f(r) = 
f(r + s). 



6 Conclusions and prospects 

This paper consists of two related parts. In the first part we have studied the representation 
theory for the quadratic algebra associated with a 2D second order quantum superintegrable 
system with degenerate potential, namely S3. We have classified the possible finite-dimensional 
representations and infinite dimensional bounded below representations, i.e., those with a lowest 
weight vector. Then we have constructed the possible Hilbert space models for these represen- 
tations, in terms of differential operators or of difference operators acting on spaces of functions 
of one complex variable. These models make it easy to find raising and lowering operators for 
the representations and to uncover relationships between the algebras and families of orthogonal 
polynomials. Here S3 has been treated as an example of a degenerate potential superintegrable 
system. The example S9 of a nondegenerate potential was treated in [39]. In 2D there are 
13 equivalence classes of superintegrable systems with nontrivial potentials: 7 nondegenerate 
and 6 degenerate. Results for all of these cases will be included in the thesis of the third author. 

In the second part of this work we have taken up the study of models of the quadratic algebras 
associated with the classical second order superintegrable systems. In each model there is only 
a single pair of canonically conjugate variables, rather than the 2 pairs in the original classical 
system. We showed, based on classical Hamilton-Jacobi theory, that such models always exist. 
Then we described a procedure to derive the one variable models for the quantum quadratic 
algebras from the models for the classical quadratic algebras. Since it is known that there is 
a 1-1 relationship between classical and quantum second order superintegrable systems (even 
though the algebras are not the same), it is not too surprising that one should be able to 
compute the quantum models from the classical models. However, we have made this explicit. 
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We applied this procedure not only to obtain the differential and difference operator models 
for system S3, but also for the generic system S9. For S9 we showed that there is a difference 
operator model associated with general Wilson polynomials, but no differential operator model. 
This construction demonstrates that the theory of general Wilson polynomials is imbedded in 
classical mechanics in a manner quite different from the usual group theory (Racah polynomial) 
approach. 

There is much more work to be done. Once the models are worked out and the corresponding 
functional Hilbert spaces are constructed, usually Hilbert spaces with kernel function, then one 
needs to find intertwining operators that map the model space to the space on which the 
quantum Schrodinger operator is defined. Also, we have demonstrated how to determine the 
classical models and show how they quantize in a unique fashion. A puzzle here is that we are 
finding classical models corresponding to non-hypergeometric type variable separation. These 
classical models typically involve elliptic functions. We do not yet understand how they can 
be quantized. They clearly do not lead to differential or ordinary difference operator quantum 
models. 

Another part of our effort is to study the structure of quadratic algebras corresponding to 3D 
nondegenerate superintegrable systems, and to find two variable models for them. This is a much 
more difficult problem than in 2D, where it led to general Wilson and Racah polynomials, among 
other models. The quadratic algebra still closes at order 6 but now there are 6 linearly second 
order symmetries, rather than 3, and they are functionally dependent, satisfying a polynomial 
relation of order 8. There are 4 commutators, instead of 1. For the models we expect to find 
multivariable extensions of Wilson polynomials, among many other constructs. 
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